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Abstract 

We consider the problem of finding four different rational squares, 
such that the product of any two plus the sum of the same two always 
gives a square. We give some historical background to this problem 
and exhibit one such quadruple. 

1 Introduction 

The genesis of this work is the following section taken from Chapter 19 in 
Volume 2 of Dickson's "History of the Theory of Numbers" p. 

Fermat treated the problem to find four numbers such that the product of any 
two increased by the sum of those two gives a square. He made use of three 
squares such that the product of any two increased by the sum of the same 
two gives a square. Stating that there is an infinitude of such sets of three 
squares, he cited 4, 3504384/rf, 2019241/c/, where d = 203401. However, he 
actually used the squares 25/9,64/9,196/9, of Diophantus V, 5, which have 
the additional property that the product of any two increased by the third gives 
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a square. Taking these three squares as three of our numbers and x as the 
fourth, we are to satisfy 

34 25 ^ 73 64 ^ 205 196 ^ 

— X H = □ — X H = □ X H = □ 

9 9 9 9 9 9 

This "triple equation" with squares as constant terms is readily solved. T.L. 

Heath found x to be the ratio of two numbers each of 21 digits. 

This section generated several questions: 

1. Is there "an infinitude of such sets of three squares"? 

2. Where did the example 4, 3504384/rf, 2019241/d come from? 

3. What is the x found by Heath and can we find a smaller value - smaller 
meaning fewer digits in its rational form? 

4. What would the solution have been with the original example of Fer- 
mat? 

5. Can we use the set of 3 squares to find a fourth square, so that the 4 
squares are a solution to Fermat's problem? 

Given a set {oi, 02, ... , 0^} such that a^aj + + a^- = □ for 1 < i < j < m, 
Dujella [2] jS| called such a set an Eulerian m-tuple, in honour of Euler 
who found the first example of a quadruple of numbers. 

To be historically accurate, however, we should note that what Euler was 
considering was sets {&i, &2, ^4} with bibj — 1 = □. This is related to the 
present problem since 

aiaj + flj + aj = {at + + 1) - 1 



2 Infinitude of Eulerian Triples of Squares 

We wish to find x, y, z such that 

2 2 I 2 I 2 2 

X y + X + y = a 

2 2 I 2 I 2 1,2 

X z + X + z = b 

2 2 I 2 I 2 2 

y z + y + z = c 
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Rather than have a comphcated analysis, we show that there are an infinite 
number of solutions with x = 2. 

Thus we look for y, z with 

Rational solutions of + 4 = □ can be parameterised as t = 4//(5 — /^), 
with / rational. Thus y = 4m/(5 — m^) and z = 4n/(5 — n^), which can be 
substituted into the third part of (1). 

We find that m and n must satisfy 

cf = rn^n^ + (m^ - 4m^ + 25)ra^ + 25m^ 

for d rational. 

Define d = Y/m and n = X/m, giving the quartic relation 

= X^ + (m^ - 4m^ + 25)X^ + 25m^ 

This quartic has an obvious rational point X = 0,Y = 5m^, and so is bira- 
tionally equivalent to an elliptic curve. Using the standard method described 
by Mordell we find the elliptic curve (with m = p/q) 

j2 = _ 2(p4 _ 4p2^2 ^ 25q^)K^ + (/ - S/g^ _ 34^4^4 _ + Q25q^)K 

(2) 

with the relation n = J/{2pqK). 

The curve has 3 points of order 2 with J = 0, which lead to z = 0, so we 
need other rational points for non-trivial solutions. Evidence suggests that 
the torsion subgroup is just isomorphic to Z2 x Z2, but this would be difficult 
to prove. If true, we need the curves to have rank greater than for solutions. 

If we experiment, we quickly find that m = 1/2 gives a curve = — 
770K'^ + U6625K which has rank 1 with generator P = (245, 2100). m = 1/2 
gives y = 8/19 and the point P gives n = 15/7 and z = 21, and it easily 
checked that {4,64/361,441} is a square Eulerian triple. 

Since the curve has rank 1 we have that integral multiples of P are also 
rational points on the curve. For example, doubling the point P leads to 
the values K = 187489/441, J = 651232/9261 and so n = 376/9093 and an 
alternative z of 13675872/413271869. There are thus an infinite number of 
Eulerian triples with x = 2,y = 8/19. 
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3 Numerical Values 



It is totally unclear from Fermat's original work where his numerical example 
involving 4 comes from, as he just states this result with no supporting 
algebra or computation. As we saw in this last section, there are much 
simpler sets which include 4. 

We now consider the system 

34 25 ^ 73 64 ^ 205 196 ^ 

— xH = □ — xH = □ x^ = □ 

9 9 9 9 9 9 



We first write 



34 25 /5 , 

— X H = — \- fx 

9 9 V3 
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which gives x = 2(17 - 15f )/{9p). 

Substituting this into the second and third equations, we find the following 
two equations must have rational solutions: 

2(288/2 - 1095/ + 1241) = 
2(882/2 - 3075/ + 3485) = /i^ 

Now, in the first equation, / = 17/15 gives e = 136/5, so, if we substitute 
e = 136/5 + g{f — 17/15), we can solve to find 

17^2 _ giQg _ 23058 
~ 15(^2 _ 576) 

and if we substitute this into the other quadratic relation we find, clearing 
denominators, that we must have a rational solution to 

= 14161/ + 731544C/3 + 28206441/ + 639486144^ + 6471280836 = G(g) 

(3) 

Since 14161 = 119^ we can attempt to complete the square, by forming 
G{g) — {ag"^ + f3g + 7)^. Simple arithmetic shows that if a = 119, /3 = 
21516/7,7 = 919177353/11662 then 

2 ^ ,2 1581221503125(13328^+22275) 



G{g) - {ag' + /3y + 7) 



136002244 
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Thus, we have a square solution when g — —22275/13328, which gives / = 
142415972261/56567733755 and finally 

_ -459818598496844787200 
^ ~ 631629004828419699201 

where both numerator and denominator have 21 digits. 

Because equation (3) has leading coefficient 119^ it can be transformed into 
an eUiptic curve. Mordell's method leads to the curve — + 20478A;^ + 
99801585A; with 

_ 3(105j - 44(163A; + 2348685) 
^~ 34(49A; + 560880) 

This curve has 3 points of order 2 with k = 0, —7999, —12483 all with j = 0. 
The first leads to a; = whilst the other two give undefined values for /. Thus 
we need points of infinite order. Cremona's mwrank program shows the rank 
to be 2 with generators (-9984, -222768) and (-8379, -114912). The first 
gives g = -543/8, / = 269/147 and x = -50176/72361, significantly simpler 
than Heath's value. 

With regard to the {4,3504384/203401,2019241/203401} triple of squares 
we can perform an identical analysis. The corresponding value of x from 
completing the square is 

_ -28448417598272924003671204878289354665765410185967616 
^ ~ 36828906078832095599985737816846193226885934523284161 

where both numerator and denominator have 53 digits. 

Attempts, as before, to find a smaller value of x lead to the elliptic curve 

v'^^u^ + 10450883424805xi^ + 26734915668323655104674200ii 

The completing the square values lead to a point of infinite order on this curve 
with u = -9390695817653070336/2019241. Investigations with APECS, 
mwrank, and SAGE were unable to find other generators. Both APECS 
and mwrank give 3 as an upper bound for the rank. The root number is — 1, 
so the parity conjecture suggests the rank is 1 or 3, but we are unable to be 
exact. 
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4 Square Quadruples 



We now consider the problem of finding {xi, X2, 0:3, X4} with all 6 combina- 
tions x'fx^ + x'^ + x^ = □ for i 7^ j. For ease of notation, we call the quadruple 
{x,y,z,w}. 

Assume x = e/f, with e and / positive and having no common factors. 
The rational solutions to (x^ + l)t'^ + x"^ — H can be parameterised hy t — 
2xm/{x^ + 1 — m^), so that 

2x171 2xn 
^ x^ — + 1 — + 1 

Assuming that m = g/h,^Q find that y'^z'^ + y"^ + z^ = D requires that 
m^n^ + (m^ - 4:m^ + x^ + 2x'^ + l)n^ + m^{x'^ + + 1) = □ 

Considering m and x as parameters, this quartic in n has an obvious solution 
when n = 0, so is equivalent to an elliptic curve. We find the following curve 
with integer coefficients, = + AK"^ + BK where 

A = -2(e"/i" + 2e^fh^ + - 4yV + h^)) 

and 

B = e'h' + 4e^ fh"" - 2e^fh\g^ + Ag^h^ - 3h^)- 
Ae'fh'ig^ + Ag^h^ - h^) + - %g%^ + lAg'h'' - %g^h^ + h!") 

withn^ J/{2ph?Km). 

We search this curve for integer points and generate n from these, giving the 
Eulerian triple {x^, y^,z^} . 

We now need to find w with (x^ + l)w'^ + = □ assuming that a; is a known 
fixed value. From before, w = 2rx/{x'^ + 1 — r^), which we substitute into 
{y^ + \)w^ + = □, to find that r must satisfy 

y V + (2a;2(y2 + 2) - 2yy + x^y^ ^ 2x'y'' + = □ 

This quartic has an obvious solution when r = 0, so it equivalent to an 
elliptic curve. We find the curve 

V^^U^ + {x\2y^ + 1) + y2)c/2 + (xV(z/' + 1) + x''y\y^ + 1))U 
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with r = V/{y{x^{y'^ + 1) + U)). 

This eUiptic curve can be easily transformed to have integer coefficients and 
then searched for integer points, some of which lead to values of w different 
from X, y, z. The final test is to compute w'^z^ + w'^ + z"^ and to test whether 
this is square. 

This methodology was coded using the simple multiple precision UBASIC 
system, and run for several hours on a PC. The code finds hundreds of 
examples with 5 of the 6 identities equal to a square, but only one example 
with all 6 square. This solution is 
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